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ABSTRACT 
Yield curves for *9Ca(d,a) reactions have been 
measured in the range 4.50 to 5.50 MeV. 2a (agp) yield 
curves were obtained from 4.88 to 5,50 MeV; both sets of 
data were taken in 20 keV energy steps. The yield curves 
were” analyzeds for“ Erieson™t luectuations, giving Mmeathe 


H2oG at 


average width of a level in the compound nucleus 
Piva tev vexcitavion, MiAwase foundetombeni 238 ikeVi 
Ansulare’distributtonsefor several *9Ca(d,p) tloa 
reactions have been measured in the range 4.50 to 5.50 MeV 
in 2 keV steps. These data were averaged over energy 
nLemveadiewot > 107.) > 20 T and > 40 1. The resulting 
energy averaged angular distributions were then compared 
with either DWBA or statistical model cross section pre- 
dictions. It was found that the two smaller averaging in- 
tervals are inadequate. Even an averaging interval of 
> 40 T is only partially successful in removing energy 


dependent cross section fluctuations, permitting comparison 


with theory and the extraction of spectroscopic factors. 
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CHAPTER el 
INTRODUCTION 
In the past decade, acs + d reactions have been in- 
vestigated by several groups (e.g. Be65, Mi65). These investi- 
gations have produced a great deal of nuclear data; in fact 
our fairly complete knowledge of most of the exit channels for 


40 : , ; 
Ca + d reactions means that a reliable estimate can be made 


of the compound nuclear contribution to various exit channels 
using the theory of Hauser and Feshbach (Ha58), (This theory 
is also known as H-F theory or the statistical model theory.) 
In other words, we can test H-F theory by comparing its pre- 
dictions with measured compound nuclear cross sections. Since 
H-F theory only predicts energy averaged cross sections, 
angular distributions of protons from *0C0a(d,p) reactions were 
obtained in 2keV energy steps from 4.5 MeV to 5.5 MeV bombar- 
ding energy. Them the measured cross sections were energy 
averaged before comparison with H-F predictions. 

H-F theory also predicts the average lifetime of 
compound nuclear states. As explained in Chapter 2, this 
quantity can be measured by considering the fluctuations in 
a yield curve; thus we have a further check on the theory. 

At the bombarding energies involved in this experi- 
ment, two different reaction mechanisms are involved: compound 
nuclear reactions and direct reactions. These will both be 
discussed in detail in the theory chapter. As mentioned above, 
compound nuclear reaction data will be compared to H-F theory; 
data from direct reactions can be compared to the predictions 


of the distorted wave Born approximation (DWBA) theory. Lee 
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et al (Le64) have in fact used the *UcaCd,p) reactions as a 


test of the DWBA. They found that the predictions of the DWBA 


for the angular distributions were very good, except for poor 


fits at back angles. It was felt that perhaps a comparison of 


energy averaged data with DWBA would be more appropriate, at 


least at low bombarding energies. 


following 


In the experiment reported in this thesis, the 
measurements were made: 

dling) Gata: a) “Oca, siyequcerc distributaons of ihe 
elastic scattering of deuterons were measured. 
Energy averaged distributions were extracted from 
data taken at 2 KeV intervals from 4.5 MeV to 5.5 
MeV. These data were used in the calculation of 


optical model parameters for use in DWBA and H-F 


calculations. 
40 41 
G2) Ca(Cd,p) Ca. Energy averaged angular 


distributions were obtained as in (1), and compared 
with DWBA and H-F predictions. Also, yield curves 
were obtained for fluctuation analysis. 

(3) +0oa(d, a) Eke Angular distributions were 
obtained at several bombarding energies. Yield 


curves were obtained for fluctuation analysis. 
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CHAPTER 2 
THEORY 
2314 ptireal,s Model. 

In both Hauser-Feshbach and DWBA calculations, we 
need to know the potentials that influence the interacting 
particles. These potentials are usually obtained from the 
optical model of the nucleus. In this model we assume that 
nuclear scattering is analogous to the scattering of light 
from a translucent sphere. 

When a wide beam of light is incident upon a trans- 
dicen t sphere. yi hemeollows neyoceuns 

(1) part of the beam is reflected; 

(2) part of the beam is absorbed and re-emitted 

with a phase change (refraction); 

(3) part of the beam is absorbed and re-emitted 

with a changed frequency; 


(4) part of the beam misses the sphere. 


themnuci#eus vs chought} ofwassey sphere of nuclear 
matter, transmitting, reflecting, refracting and absorbing 
the incident nuclear waves. Using this model we can under- 
stand many phenomena observed in nuclear reactions. for 
example, interference between the incident and refracted 
neutron waves explains the complicated broad maxima and 
minima observed in the total neutron cross section between 
.l MeV and a few MeV. The shift of these curves with in- 
creasing A corresponds to diffraction by successively larger 
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The optical model of the nucleus is then a com- 
plete “analogs to. piysicalsoptics, Ihe nucleus can be replaced 
by a region of fixed refractive index (an attractive poten- 
tial well) and a given opacity (an imaginary potential well). 

Une usta leversion ote tnesopiical potential includes 
a berm Oy account som they Spin -ompike Anieractwon.s yllis sas 
necessary toptat polarizapionsdata. ~Allso, at, the, energies 
USCG MIM ees Invest Pation es tiekr alli principle Aniubies 
collisions, particularly in the dense nuclear interior. 
Therefore, most absorption takes place near the nuclear 
surface, and our imaginary potential well must give this 


effect. The optical potential used has the form 
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(3) fws is the Wood-Saxon form factor 
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STENT ete 
1 + exp[(r-r,al/3)/al 


(ooecnded sane radius and diffuseness 
parameters, respectively) 
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(4) V is the strength of the refractive part 

of the potential. 

(5) W is the strength of the imaginary part of 

the potential. W accounts for absorption of 

particles from the incident beam, i.e. attenuation 

of the beam. Note that W is multiplied by the 

derivative of the Wood-Saxon form factor; this 

gives the desired peaking of W at the nuclear 

SuUrLaCe. 

(6) Vso is the strength of a spin-orbit inter- 

action of the Thomas type, which is also peaked 

at tne ssurtace., 

(7) M,_ is the pion mass. 

ihewsnapes (Om Vea Wi and VSO, (aS, well as the real 
part of the radial wavefunction for %£=0 incident deuteron 


anew shown on fig. 2... Lhis as taken from Bot. 
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Figure 2.1 Typical shapes of real, imaginary and 

Spin-orbit potentials, labelled V, W and Vso respectively. 
The curve labelled "incident deuteron" shows the real part 
of the radial wavefunction for an & = 0 deuteron incident 


on the shown potentials. 
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2.2 DWBA 

As mentioned in the introduction, at the bombarding 
energies used in this experiment two types of reactions are 
significant--direct reactions and compound nuclear reactions. 
In this section we will consider direct reactions. 

inva direct reaction, the transition from the in— 
cident channel (initial state) to the reaction channel (final 
state) is considered to take place in one step without an 
intermediate step. The reaction typically takes place in 
about 10722 sec., the time it would take the incident deuteron 
to pass the diameter of the target nucleus. A special case 
OF the direct reaction as the stripping reaction. In a 
stripping reaction, the target nucleus captures part of the 
ineident™ particle, but not all ofr it.” The remainder of the 
particle leaves the vicinity of the target without causing 
further major effects, Thus a direct (d,p) reaction is a 
stripping reaction. 

Consider a stripping reaction A(a,b)B, where 

A is the target nucleus 

a is a deuteron 

b iS a neutron or proton 

Be isitthe residual nucleus, 

The calculation of the transition amplitude between 
the. initial state (A,a) and the final state (b,B) is performed 
tor the tirst order only in the interaction potential between 
the incident particle and the target nucleus. Then, using 


perturbation theory, we obtain the following expression for 
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the transition amplitude; 
ieee Wea ae eles bid vole es Cae 
Where $¢¢ is the asymptotic wavefunction of (b,B) 
¥; 1s the initial total wavefunction of the system 
V* is the perturbing potential. 
(Recall that the cross section is proportional to |T|?) 
The calculation of T can be greatly simplified by 
assuming ¥; is the product of an incident plane wave and the 
internal wavefunctions of A and a. This is known as the Born 
Approximation. Of course, this approximation ignores the 
distortion of the incident plane wave by the repulsive 
coulomb force, and the nuclear interaction. If we calculate 
T using the more appropriate distorted waves, we are making 
the Distorted Wave Born Approximation or DWBA. These distorted 
waves are simply the wavefunctions found from solving the 
Schrodinger equation 
v2) Cp) yen CE > ylr)]y(e) = 0 g.e+---22+--2-3 
where u is the reduced mass, U(r) is the optical 
potential, described in the last section. 
Since U(r) is obtained from the analysis of elastic scattering 
measurements, and since the interaction potential V~ in equa- 
tion 2.2 is treated as a perturbation, we are implicitly 
making the assumption that in our reaction elastic scattering 
predominates, and the other channels make only small con- 
tributions. 
Now, for a stripping reaction of the type A(d,p)B, 
the total Hamiltonian of the system can be written as 
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where TdA is the operator for the relative kinetic 
energy of the deuteron-target system 
VpA is the potential energy of the proton- 
Lange vesySten,eecier 
The Hamiltonian for the final state, He, can be 
written 
nie ae WW Ge Abels) oe havi, A Wie AAA eG 5g oon es 
where VpB is the optical model potential determined from the 
elastic scattering of protons on B. From equations 2.4 and 
2.9 we see 
Ng Dig pee DE. 
where V’is the perturbing potential in 2.1. 
We now make the assumption that 
VpA = VpB 
and we have 
Nie oeea Dlvaelenwtels 6 sls titekeushe ore sregener cis /0 
(This derivation for V- follows that of En66)). 
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T = <¢¢|Vpn| ¥y> 


Ftyy* (x, 75) <B| Vpn |Arxg "ka stg drparg 
where x, and Xp are the incoming deuteron and outgoing proton 
distorted waves, respectively, |A> contains the deuteron 
internal wavefunction and the target nucleus wavefunction, 
<B| contains the residual nucleus wavefunction. We can 
simplify the integration by assuming that B= Atn. Then when 
we calculate the overlap of the internal wavefunctions of A 


and B, only the wavefunction of the captured particle remains. 
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nV onbd 
Here $n is the wavefunction of the neutron bound to the 
targer nucleus A, and ¢q as the internal waverunctson of the 
deuteron. 

A further approximation is now made: the range of 
VG toese GtG sZero, so Enemproduen Von?d Can) Der rep aceam yy. 
deta function multiplied by a constant, “This implies thar 
the proton appears exactly where the deuteron disappears. 
This is unsatisfactory because the deuteron is a large, 
diffuse particle. However, we can simulate the effect of 
the finite size of the deuteron by multiplying the zero range 
form factor by a constant correction factor (Buby), (Peo4). 
This approach is known as the local energy approximation, 
and the correction factor is called the finite range 
COLVreetTUCh. 

Nuclear forces are thought to be non local-—-that 
ae the Particle interacts with the potential over all space, 
NOG Only at the center of its wave packet. Thus we should 
replace w(r) Ulr) in the Schrodinger equation (2.3) by 

fu(e"-)u(r,r-)d*r-. However; to facilitate the calculation 

we simply multiply the wavefunctions of the local approxima- 
tiLGneby va, constant sdependingponsthe sreal ywelledeptheandgthe 
range of the monlocality CHj 65). 

Thus, sity oun direct intenactionyjcalculations ewe 
have made the following assumptions: 


(1) That the Distorted Wave Born Approximation 
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is valid for our experiment; 

(2) That V = Vpn (see equation 2.6); 

(3) That the internal wavefunction of the residual 

nucleus is the wavefunction of a neutron bound to 

the target nucleus. The target nucleus is treated 
as an inert core. 

(4) That the local energy approximation is valid; 

(5) Non-local effects can be approximated by mul- 

tiplying our local wavefunctions by a constant. 

All of these assumptions are thought to be reason- 
able, and in fact are the standard assumptions made for the 
analysis of stripping reactions. 

It is seen in assumption (3) that we are coupling 
the. stripped neutron to an inert 4Oca core. The probability 
that for a given state the configuration of the nucleus is 
that of a single nucleon orbiting about an inert core is the 
Spectroscopic) factor Tor) that state. Therefore our calcu- 
lated cross sections correspond to a spectroscopic factor of 
1. We can measure the spectroscopic factor S of a state by 


finding the ratio 


g = (do/d2) experimental 
(do/d2) DWBA 


The shape of an angular distribution calculated 
using the DWBA is strongly dependent on £, the number of 
units of orbital angular momentum exchanged between the in- 
cident particle and the target nucleus. Thus a comparison 


of experimental and theoretical angular distributions can 
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give us &. Furthermore, if the target nucleus has spin and 
parity 0t, once 2 is determined we know 


Je = Mt. 1/2 
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2.3 Compound Nuclear Model 


As mentioned before, we are concerned with two 
types of reactions: direct reactions and compound nuclear 
reactions... Direct reactions were discussed in section 2.2; 
this seetionsdeals with compound. nuclear,;reactione. In this 
type of reaction, the incident particle is completely absorbed 
by the target nucleus, sharing its energy with some or all of 
the target nucleons. Eventually, enough of this energy will 
be acquired by one nucleon or group of nucleons to allow it 


to escape the system. The time interval between absorption 


of the incident particle and emission of the outgoing particle 


is relatively large--typically 107-1© seconds when a neutron 
is absorbed. This means that a fairly long-lived entity has 
been formed; we call it the compound nucleus. 

Reasonable predictions for compound nuclear cross 
sections can be made by making the following arguments (see 
Vo69 and references therein). 

We consider the reaction 

a + CN + a” 
where a labels the target nucleus, the incident particle, 
and their state of excitation. (% is called the 
entrance channel.) 
CN represents the compound nucleus. 
a’ labels the residual nucleus, the outgoing par- 
ticle, .and.thein state of excitation. %° 15 one 
of the exit channels. 


We assume that 
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14 
(1) The cross section can be factored into a 
part representing cross section for formation of the 
compound nucleus in the entrance channel, and a part 
representing the probability that the compound nucleus 
will decay into the appropriate exit channel. 
(2) Total angular momentum and parity are 
conserved, 
(3) The optical model can be used to calculate 
the compound nucleus formation cross sections. 
(4) The following reciprocity relation can be 
wsec: 
Ol r= seo Te 
where k is the wave number for the relative 
motion of the pair a 
G.. 12 LS the cross section for the reaction a>a' 
Tl, 2s the cross section for the reaction a!+a 
(5) The random phase approximation can be made. 
Essentially this means that the compound nuclear states 
are so complicated that the matrix elements connecting 
a group of them to any other given state are random 
complex numbers. The matrix elements vary rapidly 
with energy, both the real and imaginary parts assuming 
positive and negative values with a mean value of zero. 
Now, interference between different partial waves is 
manifested in the theoretical cross section by terms 
involving the product or difterent matrix elements. 


Therefore, if we perform an energy average these inter- 
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ference terms will cancel because of the statistical 
fluctuations of the matrix elements about zero. This 
cancellation of interference gives the characteristic 
Symmetry of the compound nuclear cross section about 
90°, as if we were dealing with isolated resonances 
of definite parity in the compound nucleus. 

Because Of assumption (6), the model desenibed 
above 15 called thesstacistigal gmodelgarine model is also 
called the Hauser-Feshbach model, after the authors who 
first added conservation of total angular momentum and 
parity to the “theory (ha5s2)y. 

With the above assumptions, we obtain the follow- 


ing expressions for compound nuclear cross sections: 
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where Gabe is the energy averaged total cross section 
I and i are the intrinsic spins of the pair of 


particles a 


S usyihe channel spin: 
Ss 1p 
Q is the orbital angular momentum of the pair 
J is the total angular momentum: 
ee) ee 
Il is the total parity (the product o: (—)* 
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a, given by the optical model phase shift éa2: 
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Myrcodes 1 -« |e2i6eh|2 
The sum over unprimed quantities is a sum over 
entrance channels, and the primed quantities are summed over 
exit channels. The double primed quantities are summed 


over all possible decay channels. 


For differential cross sections 
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where the Z's are angular momentum coupling 
coefficients defined by Biedenham, Blatt & Rose (Bi52). 

In both 2.7 and 2.8 we have sums over all possible 
exit channels in the denominators. In this investigation, 
decays leaving the residual nuclei in highly excited states 


OED iB Since there is little 


were permitted (e.g. 11 MeV for 
spectroscopic information available for such exit channels, 
we are forced to replace the sum over discrete channels by 
an integration, which requires an estimate of the density 
of levels as a function of energy. This will be discussed 
in detail in section 4.5. 
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average our data over a large enough energy interval for 
the random phase approximation to hold. If this is not done, 
we cannot expect either symmetry of our compound nuclear 


eross sections about 90°, or quantitative agreement between 


measured and predicted cross sections. 
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From quantum mechanics we recall that 

(1) The probability of an event in an experiment 
is given by the absolute value of a complex number 
¢ which is called the probability amplitude: 


P's. probability 


$ 
o 


probability amplitude 
lo |? 


(2) When an event can occur in several alternative 


ways, the probability amplitude for the event is 
the sum of the probability amplitudes for each way 
considered separately. There is interference: 
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Thus the overall probability is the probability 
that the event occurs by mechanism 1 plus the probability 
that the event occurs by mechanism 2, Are a term indicating 
interference between the two mechanisms. (This treatment is 
from reference Fe65). 

In this experiment, we are measuring cross sections 
of states that can be formed either by direct or compound 


nuclear reactions. In this context equation 2.9 can be 


written: 
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where On is the total cross section 


Scy is the compound nuclear cross section 
On iS) tne direct cross section. 

In the most general case, then, part of the total 
cross section comes from the interference between the direct 
and compound nuclear mechanisms. However, if we obtain 
energy averaged cross sections so that the random phase 
approximation is valid, the phase differences between the 
direct and CN amplitudes will be random numbers, fluctuating 
rapidly with energy about a mean of zero. Thus, the inter- 


ference terms cancel and we are left with 


on = Soy + 9p 

where o is an energy averaged cross section. 

It should be noted that it is still debatable 
as to whether or not interference between direct and CN 
amplitudes has ever been observed. Many investigators 
Simply assume that direct and compound nuclear reactions are 
incoherent, i.e. that interference is negligible. This 
assumption probably introduces a much smaller error than 


the errors inherent in DWBA and Hauser=-Feshbach calculations. 
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As explained in section 2.3, if we measure a 
cross-section as a function of energy, we will observe 
fluctuations above and below the mean cross-section unless 
the random phase approximation holds. Such a yield curve 
is made up of overlapping compound nuclear resonances 
which may interfere constructively or destructively, giving 
pEse to the fluctuarirone: 

Ericson noted that these fluctuations could be ex- 
ploited to give information about compound nuclear properties 
(ER63). If the overlapping resonances have an average width 
T, a particular resonance will influence the yield curve over 
an energy range of approximately I. Now compare the cross- 
section at two energies, E and E + e; 

(1) if e<<I, then at both energies the same 

resonances contribute to the cross-section, and 

the two cross-sections will be quite similar; 

(2) if e>>T, different resonances predominate 

at the two energies, and the cross-sections will be 

lnecopre laced. 

Thus, by measuring the correlation of the cross- 
SeCclLicone Ss vaarUnNctLOnN Of te, We cam estimate f.. = lnevcormelaq 
tion 1s measured using the autocorrelation function: 

aoe qa Ce) =O io CEte) = cil 

Co) 
where < ... > 1S an enerpy average 
E is the energy 


o is the average cross section. 
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Once [ has been measured, we can obtain the average 
Compound nuclearndafet imerasingethel uncertainty. pninciples 

Pr oth 

where t is the average CN lifetime. 

qt can be predicted by a Hauser-Feshbach calculation 
and compared to the measured value. 

Often a fluctuation analysis is done on a state 
that has a significant direct reaction cross section. Thus, 
the fluctuations due to CN reactions are superimposed on a 
direct reaction cross-section. Ericson derived a form of 
Thesautocorre lations function, for suchria. situation: 


ie) ee gi nt ee 
N DE BAe 


where N is the fluctuation damping coefficient (see 
below) 
On ie the Yatiovor direct, to totalveross— 
section. | 
Since we are using an unpolarized beam, contribu- 
tions from different spin projections are incoherent. If 
there are N equivalent inchoherent contributions, the fluc- 
tuations in the yield curve will be reduced by 1/ YN. Near 


90°, an approximate value for N is 


N = & [C2T+1) (2941) 027141) (23141)) 
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Unfortunately, the finite range of data errors (see below) 
made it impossible to extract meaningful values for O5 from 
this experiment. To obtain ie the yield curves could have 
been taken over a much larger energy range; time limitations 
made this unfeasible, 

The main error in a fluctuation analysis comes 
from the finite energy span of experimental yield curves. 
In other words, if more CN resonances enter into our ex- 
perimental yield curve (i.e. if we extend the measurement 
over a larger energy range) we can determine I more 
accurately. Expressions estimating the errors introduced 


by the finite range of data have been derived (Pa65). 
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CBAPTER 3 
EXPERIMENTAL 

3.1 Scattering chamber, detectors and associated equipment. 

The incident deuteron beam for this experiment 
was produced by the University of Alberta CN type Van de 
Graaff accelerator. The measurements were taken in a 16" 
scattering chamber designed by Dr. G. Roy, Five silacon 
surface barrier detectors were used, simultaneously collect- 
ing data at different angles. Circular collimators were 
placed in front of each detector. They subtended a solid 
angle of about 2 x 1077+ steradians. The relative solid 
angles of the collimators were measured by placing each de- 
tector at a given angle in turn and observing protons 
elastically scattered from gold. To decrease noise and im- 
prove resolution in the detectors, they were cooled with an 
ethanol-dry ice mixture. For the angular distributions, four 
movable detectors were mounted 10° apart at a distance of 8" 
from the center of the chamber. A fifth detector held sta- 
tionary at 90° was used as a monitor. For the yield curves. 
five stationary detectors were mounted 8" from the center. 
The voltage pulses from each detector were processed by a 
pre-amplifier, an amplifier and an analog to digital con- 
verter (ADC) in turn. That is, there were five independent 
systems, consisting of detector--pre-amp--amp--ADC. Each 
ADC was connected to a clock for dead time measurements. 
The output of the five ADC's was stored an the data acqui-— 
Sition computer (a Honeywell 615). “At the end) of a run, 


the data were transferred to the data analysis computer 


(an SDS 920) for writing the data on magnetic tape, summing 
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23 
Peaks, plotting (spectra ete. “Gsee fic. 3c.1)." “As soon as 
the data were transferred to the 920, the 516 was free and 
a new run was begun. This saved many hours of machine time, 
Since data collection, storage and analysis could be carried 
out - the same time rather than sequentially. 

The beam was collimated by two adjustable beam 
defining slits and an anti-scatter slit, (Ro69a). The beam 
defining slits were adjusted to minimize background counts. 
After passing through the target and out of the scattering 
chamber, the beam was collected in a Faraday cup. A secon- 
dary electron suppressor biased to - 1000V was mounted 


Slightihy in. trout. of the Faraday cup. 
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3.2 Target biasger. 

As mentioned in the theory section of this thesis, 
experimental compound nuclear cross sections must be energy 
averaged before comparison with Hauser-Feshbach predictions. 
In other words, we must obtain cross sections at a number 
of bombarding energies in a large enough energy range 
that the random phase approximation holds, and then average 
the cross sections. There are several ways in which data 
can be obtained for a range of bombarding energies: 

(1) If a very thick target is used, the beam 

energy is degraded over an appreciable energy 

range. This method is used very frequently in 
gamma ray spectroscopy--e.g. see Pi/70. However, 
outgoing particles are also degraded over a large 
energy range if they can get out of the target at 
alijeso this ,method.isSenot, acceptable,~for,hagh 
resolution charged particle work. 

(2) The terminal voltage of the Van de Graaff 

can be changed. In this method we measure cross 

sections at one energy, Slightly change the bom- 

barding energy, measure the cross sections again, 
and so on. After the run the cross sections are 
averaged manually. This is obviously very tedious 
and consumes large amounts of machine time. 

(3) -A.voltagescan, be applied, te the tansetate 

slow down the bombarding particles. If the out- 

going particles have the same charge as the incoming 


particles, they will arrive at the detectors with 
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almost exactly the same energy for all values of 
the bias voltage applied to the target. For 
example, by applying a bias voltage continuously 
varying between 0 and 50 kv during a run, we can 
simulate the effect of a thick target, but retain 
good resolution by actually using a thin target. 

In this investigation, a combination of (2) and 
(3) was used for the (d,p) work. During a run, the target 
bias was stepped between 0 and + 50 kv in two kv steps. 

After an angular distribution was completed, another was 
taken 50 keV lower in bombarding energy. Essentially, we 
have taken *Uca(d,p) angular distributions every 2 keV from 
4.50 to 5.50 MeV bombarding energy. We have automatically 
averaged cross sections over 50 keV intervals, and manually 
averaged the 50 keV averaged cross sections over larger 
intervals. Since deuterons are singly charged and alpha 
particles are doubly charged, this technique cannot be used 
for high resolution (d,a) work. 

The target biasing system consisted of a 0 to 60 KV 
de power supply connected to the target. An automatic con- 
troller for the system was designed by Dr. G. Roy and Mr. L. 
Holm (Ro71). In the automatic mode, the system functioned in 
the following way: 

The Van de Graaff terminal voltage was set for a 
particular deuteron bombarding energy, and the bias on the 
target was set to zero. The run was begun, and the amount 
of beam passing through the target was monitored by measuring 


the amount of charge collected at the Faraday cup. After a 


yet ot bobiqas gst 
iincnini hes Pyar ais. 6 Ne NSee Ae pe? he 
aed be /, ning 8 AER vt 0?) aa ‘meawred gatvakv “AP a a 
niptex tod ,fegrsy Aolds 5 to: sostte oft otslumte . 7 | 
~Isg¢es Mist p Aotaw ylisutss yd roituloest boos 
bis (%) to Nokteatdmos s ,noktsgitesvnt sandy ol ad 
ssgust oft nor 5 gaiwad .Atow (a,b) eit xot beev asw (6) r 7 


.eqats vA owt mi va 02) + Bre 0 raswied bsggeTe abw ake 
auw verftots ,betelqmon asw noitudiatetbh asivyae oe goeRA 


$a. Vilsites2e!7! \yprsis garbiadmod ni yswol Vet 0c nsaey 


mowt Vet S yisve enmottudizteib selugns (a, b3s08" nsAst eved @ 
vilécttsmetus Sven sw WWeysns gnibsecriod VeN Beye oF o2.e : 
yifsuism bre ,elevistalt Vsx Od tove snottises fe2onro begsxsvs 
tsqisi “svo eaoktose e20ts *boxsiavs Vor Uc sft bogstevs | 
giqis big beg tad> ylgnse 5x6 Boorssusb sonke .aisvestnk : 
beau ed tonaso supindioat ebdt .baagrsdo’ ySdvob sis eslobtamg : 
: 


<ANow (a,b) noisuvlosset dgked vot. 

YA 02 oF DO f to bstarxnop méteys agnitesid toques SAT 
=nOo sitamotus oA “;tesi6s sit oF betosnnod yiqque tswog OB a 
jf (2M Sune woR 1D sad yd benmaresb asw meteye oft to? wetlogd 
fi bSnokTsnut medeye eAt ,sbom oigsemotus sit al .(1%oA) mick 
:¥6W gniwollot sat - 

5 101 tee usw sqetioy Lenimtet Vissi) sb msV off : 

oft Ao geld otis Sas .Ygions’ gntbysdmod nodsiush wslustzasq: | 

SawoRs ot bas pouged aewinis SAT LoNBS' Os 352 cow tome 
acento bexotinom eaw tegrar ‘ie tao uaieosq masd 20 | 
BRT que vsberet oy ap batintest nieete 10 sayoms oxf3 


2] 
specific amount of beam had passed through the target, 


the control system temporarily stopped the run, incremented 
the target bias by + 2 kV, and restarted the run. The bom- 
barding energy was then 2 keV lower than previously. When 
+ 50 kV bias was reached by this process, two runs at that 
voltage were taken; then the bias voltage was automati- 
Cally decreased by 2 kV per step. The run was continued 
until 0 kV bias was reached again; data were taken at 0 kV 
and then the run was terminated. Thus, we made two measure- 
ments at every target bias voltage between 0 and + 50 KV 
inclusive. The cross sections measured from the resulting 
spectra were therefore already averaged over a 50 keV 


energy range. 
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3.0 Tlapgets 

Targets were made by evaporating natural calcium 
onto backings of thin carbon on formvar. (See table 3.1 for 
the isotopic abundances of natural calcium.) The carbon 
layer was necessary to dissipate the heat generated by 
the beam; it was kept very thin in order to minimize its 
effect on energy resolution. The formvar layer was necessary 
for mechanical strength. The calcium was evaporated onto 
the formvar rather than the carbon, since it has been found 
that the calcium apparently reacts with the carbon and 
weakens the target (Bo69). 

About 30 ug/cm* of calcium was deposited on the 


backings. This gives an energy loss of <5 keV in the target. 
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3.4 Absolute cross sections. 
experiment, we must obtain a factor that will convert 
relative cross sections (e.g. the number of events observ 
ina particular reaction, suitably normalized) to absolut 
cross sections (e.g. the cross section for a particular 
reaction in millibarns per steradian). This was done in 
Slightly different ways for the Caine experiment and 
4+O0Ca(d,p) experiment. 
(a) *ScaCd,a). The elastic Scattering of 2.0 
ProLons from +Oca was Wea suredectero Ug. 0 geamU. cnc oe 
degrees. This reaction proceeds almost exclusively 
Dy the Rutherford scattering mechanism. Since the 
Rutherford differential cross section is a function 
of the bombarding particle's charge and energy, the 
Ciarege. Of ulle tapes nucleus, and the angle oF 
observation only, it can be easily calculated. 
Then the number of counts observed at each angle’ can 


be corrected for electronic dead time and collimator 


millibarns/steradian 


i nd the a or a epiven 
size, and ETO hormalized count f E 


amount of beam passing through the target can be de- 
termined. This factor can then be used to convert 


the number of counts to absolute cross sections for 
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In order to measure absolute cross sections in an 


ed 


S 


the 


MeV 


other reactions at other energies, providing the same. 


normalizations are made for collimator sizes, dead 


times and the amount of beam. As a check, during the 


(d,a) angular distributions occasional runs were made 


at 5.00 MeV, and the 4Oca(d,d) elastic peak was summed. 
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The absolute (d,d) cross section was then compared 
to the values reported by Leighton (Le68). These 
check runs were always consistent with Leighton's 
results, within experimental errors. The (d,a) data 
were obtained in a series of five runs, all using the 
Sameetarcet. | The *Uca(p,p) measurement was Made ai ter 
each of these runs. No evidence of target deterioration 
was found. 

Cb) pCa oe These data were obtained in eight 
runs, using four different targets. The absolute 
cross sections were determined in much the same way 
as in (a). The only difference was that for check 
runs, data was taken at 5.00 MeV with no target bias, 
and the foCatd.p) 41¢q9-0 peak was summed, and compared 
to the work of Leighton (Le68). This peak was not 
observed in the d,a® work, since detectors with a thin 
depletion zone were used for that work: i.e. the high 
energy protons were not stopped in the active region 
of the detectors. However, this peak is actually 
preferable for a consistency check since it is ina 
region of no background and has a fairly high cross 


section at forward angles. 
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CHAPTER 4 
RESULTS AND DISCUSSION 
ele fr iuctuat ton, Analysis, 

Using the experimental configuration described in 
section 3.1, yield curves were taken for the *O05(4,p) and 
*OCa(d, a) reactions. The (d,a) yield curve was measured 
from 4.50 to 5.50 MeV, while the (d,p) measurements were 
in the range 4.88 to 5.50 MeV. (The diserepancy in the two 
ranges was unplanned. It occurred when the last +a target 
broke near the end of the time allotted for the measurement. ) 
Both yield curves were taken in 20 keV steps. Measurements 
were made at 70, 100, 120, 150 and 170°. Representative 
results are shown in figure 4.1. 

Using these data, the autocorrelation function 
was calculated for each excitation curve. Recall from 
Chapter 2, part (5), that the autocorrelation function C(e) 


is given by 
Cle) spo CE) = cilioh Bre So i 
Co) 
where <...> 1S an energy average 
E is the energy 
So is the average cross section. 
Generally, the Ericson fluctuations will be super- 
imposed upon coarser structure in the yield curve. Thus, 3&3 
is a function -of energy. In this study, o was calculated 
using the "sliding average" (SA) method of Pappalardo (Pa64). 
In the SA method, the experimental yield curve is broken up 


into energy intervals 4. o is calculated for each interval 


and the entire set of values of o as a function of energy 
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Figure 4.1 


A) *06aCd,o,) 98x yield curves. Lab angle = 100°. 
B) *Ocacd,p)ttea yield curves. Lab angle = 100°. 


In both cases, the statistical error is less than 


the size of the data points. 
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Figure 4,2 

Autocorrelation funetion with © = 0 
ase a function of A, the averaging interval, The 
Caleudation is for the mK ground state yield curve 
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Figure 4.3 
Autocorrelation function with “A = 350 keV as 
a function of «. The caleulatidon is for the same data 
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is obtained by interpolating between these values. 6 is 

a function of A, thus the autocorrelation function is also 

a function of A. Figure 4.2 shows the behavior of C(e=0,A) 
with change in A. A value of A is selected on the "plateau" 
where CCe=0,A) does not change with A. C(¢) is then cal- 
culated varying «. Typical results are shown in figure 4.3. 
The width at half maximum of this curve is Tf, the average 
width of a compound nuclear level. 

r was calculated by this method for each observed 
state in the residual nuclei, at each angle. After taking 
an average of the r's obtained, the average width of a 
level in 42o, at *15 MeV excitation was found to be 18+5 keV. 
(The error quoted is the finite range of data error - see 
Chapter 2, section 6). 

It is somewhat disturbing to note that the value 
for r is essentially equal to AE, our energy increment in the 
yield curve. Although some yield curves have been analyzed 
with 4E51 (e,g. Bi70), it seems preferable to have Aecl. 
iG epee suggested that measurements with Abel tend ~to 
erroneously give values of r=AE, (Da69). Because of this, 
it is felt that the quoted value for Tf must be viewed with 
some suspicion. However, an upper limit for rf has definitely 
been established: [< 23 keV. This information is very 
useful in determining the energy range over which measured 
cross sections must be averaged before comparison with the 
statistical model. This will be discussed further in sec- 
BLOM MOreciromehaL lel. 
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Hauser-Feshbach calculation. This involves evaluating 


the following integral: 


Toe f D NGY, Esse) W, (E,-€ de 
Wo Ey MP ree ns ear ie. Miata 
where v labels the type of outgoing particle 
E, = excitation of compound nucleus 
ef =Senergy Or outgoinp particle 
W, (EB, )= Vevery density of Compound = nucleus at 
Ey Mev “excitataon 


Wy (Ee) =" level density of residual” nucleus at 


Ee MeVwexcrtation 
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td 
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cross section of inverse reaction, i.e. 
capture of the appropriate particle 
incident at e MeV on the target excited 
EOen > eicy- 

T was calculated using the program HAUSERT. All 
of the quantities in 4.1 are well known except W,, the level 
density of the compound nucleus 42sc, Level densities are 


passed to HAUSER in the following manner. Let N(E,) be the 


number of levels below an excitation ‘energy EO.) “Thenva splice 
Ole log aNCE) vs. / Bx is close to a straight line, the 
deviations becoming less and less as E, increases. The 
slope and intercept of the line determine the level density; 
these quantities are given to HAUSER. 

Figure 4.4 is a plot of log N(E,) vs. v Ey for 
42sq, Fifteen levels are known in this nucleus to an 


excitation of 3.36 MeV. Clearly, this is not enough to 


+Obtained from W. R. Smith of the Oak Ridge National 
Laboratory and modified by N. E. Davison of this laboratory 
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define a unique line. Limiting cases can be estimated 
however, and lines 1 and 2 in figure 4.13 correspond to 
T's of 15 and 40 keV respectively. To summarize, we find 
theoretically 15 << 40; experimentally [<23. 

Tt should be noted at this time= that) WwW, does nor 
enter the Hauser=Feshbach calculation of the cross sections 
for formation of particular states in residual nuclei. 
Therefore our inadequate knowledge of this parameter only 
affects our predictions for [; absolute cross sections can 


be predicted much more accurately. 
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Figure 4.4 

Level density for 4250 (the compound nucleus in 
this study). The two curves (1 and 2) through the data 
are of the form log N=a+t+b YE. Both give acceptable fits 
to the data; this illustrates the difficulty in determining 


the level density of +2Sc from the existing data. 


Wye STU) i 
at aueloun bavogmas otis) oe** 262 wrtansh Isyveud 
sisbeont Myuows (S bma 1)/ sevqus owt dT . Cvbute abel 


atit sidsfqeoos eviy dtod av d+ 5 = gol miot sit to om Z a 
4 


aninimisteb ni yskuoktti6 anit astetteulit aidt zs7sb edt oF 
«Stab gattetxse oft. Taek? oe*" to ytbensb fsvet ott a 


7 


30 if 


42 
20 a 


NCE) 


YE in (Mev)1/2 


39 


4.2 Selection of Averaging Intervals. 

Recalle.rom €nepter 2separt Vo) tiat., ror 
comparison of measured compound nuclear cross sections 
with Hauser-Feshbach predictions, the data must be energy 
averaged over a large enough energy interval for the 
random phase approximation to hold. This will also 
eliminate any direct-compound nuclear interference (Chapter 
2,7 part 4). It is*often assumed= that an interval Mor-=7UET 
TomeUrricient Lor this Cerra ot) 0) eine titer layemr teh — 


Vals OL 250 KeV C2 VOT) .o00™ kev (<2 20r)"and 41000 kev 


(2 40T) were used for purposes of comparison. 
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4.3 Optical Model. 

Optical model penamerere for the +9ca+d system 
were obtained from the measured 4 0ca(d,d) +0ca elastic 
scattering using the computer code SNOOPY 2*. This program 
will search over up to six of the optical model parameters 
to give a best fit to the experimental data. 

From the nature of the optical model, we expect 
the parameters to vary slowly and smoothly with bombarding 
energy. It has been found that the optical model para- 
meters that give the most satisfactory results for *8ca(d,p) 
*3ca DWBA calculations are energy averaged parameters. (Bo69). 
That is, the elastic scattering is measured at several 
energies, several sets of optical model parameters are 
extracted, and these parameters are averaged.. These 
average parameters are then used for DWBA calculations at 
all energies. In Bo69, this gave superior results to DWBA 
calculations using parameters obtained from elastic scatter- 
ing measurements only at the energy in question. This 
Suggests that the measurement of elastic scattering cross 
sections at one energy is insufficient to obtain good 
optical model parameters, even when they are to be used 
for DWBA calculations for reactions occurring at that same 
energy. For this reason, the parameters were obtained from 
elastic scattering data that had been averaged over 1 MeV. 
(The data were taken from 4.50 to 5.50 MeV in 2 keV steps). 
These average cross sections are presented in figure 4.5. 


SNOOPY varies the optical model parameters to 


*Obtained from W. Haeberli, University of Wisconsin, Madison, 
Wisconsin. 
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Figure 4.5 

s eaicd.d) elastic scattering angular distribution. 
Data is averaged over 1 MeV with the mean bombarding 
energy = 5.00 MeV. The theoretical curve is generated from 
the optical model. The parameters used in the calculation 
are given in table 4.2. Statistical errors are less than the 
Suzeroie wehneepounts. 

A) differential cross section vs. angle 

B) ratio differential cross section/differential 


Rutherford cross section vs. angle. 
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42 
give the best fit of optical model theoretical cross sec- 


tions to the experimental data. The "goodness of fit" is 
measured by y?: 


x2 = : (Scalc ti“ texp Pi) |? 
= Negra Shy 


where the sum is over all experimental angles 
Oo (¢.) 1s the optical model differential cross 
Cale + 
section at angle dae 
Fern 6? i) is the corresponding measured cross section. 


oper at) fs Une er VOre ny CMC eMeasUved sCross™ = cGUlOl. 


inchs Study, (¢5) is the differential elastic 


Cexp 
cross section, minus the Hauser-Feshbach contribution. 
SNOOPY generates sets of optical model parameters by mini- 
mizing x* while searching over the optical model parameters. 
Since the search is carried out in a 9-dimensional parameter 
Space, it is very easy for the program to become "trapped" 
in local minima. To avoid this, the sequence in. which 
parameters are searched, and the limits of the search must 
be chosen very carefully. In this study, the method 
suggested by Bogaards (Bo69) was used: 

(1) Vso, rso and aso were fixed at reasonable 

values; 

(2) V and a; were stepped through a reasonable 

range of values; 

(3) r_, a, wand pr, were searched simultaneously 


O 

For cach) Valuecot \@andeas., “theystantainp varues 
for the searches were the parameters given by 
Schwandt and Haeberli (Sc69), and are summarized 


in Table 4.1. 
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“Table 4.1 


Initial values of optical model parameters for search 


Po a Wo ry Vso rso aso 
fy OSer mee 0S oe im es. 5gMeve A966, meeCeMeVvetorourmoend .6. fm 
V was stepped from 80 MeV to 125 MeV in 5 MeV 
increments, while abVvarieds Fromeo L208 tor. soe ime ins 02/05 im 
increments. (Schwandt and Haeberli found 112 MeV and 0.52 fm 


Tor Veand®a. ) 


Table 4,2 gives the results of this procedure: 


Taber tre 
Optical model parameters generated from 1 MeV averaged 


elastic Beatcening data. 
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Figure 4.5 shows the cross sections calculated 
from the optical model using these parameters and the 1 MeV 
averaged elastic scattering cross sections. In this study, 
calculations were done with two sets of optical model para- 
meters: the above set, and those of Schwandt and Haeberli. 


These are compared in table 4.3. 


Table k..3 
Optical model parameters used for calculations. 
Set I--generated from 1 MeV averaged elastic scattering data. 
Set II--from Schwandt and Haeberli (Sc69). 
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44 CaCd,a) K reactions. 


IGS): or 


K angular distributions were obtained 
Ste. Us teon oe OU, so.2o and ay. sU Mey.  oInce thin targers 
were used, and no averaging was done, these distributions 
are not expected to agree well with theoretical predictions. 
The experimental angular distributions are presented in 


figures 4.6 to 4.10. No theoretical calculations have been 
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Figure 4.6 


*OCa(d,a) 38xe5 angular distributions. 


Figure 4.7 


7 OCa(d.a) 380.13 angular distributions. 


Figure 4.8 


*OCa(d,a) DR oS angular distributions. 


Figure 4,9 


soGa(as a) ko angular distributions. 
Figure 4,10 


*Oca(d,a) 38x2-41 angular distributions. 


In all cases, the statistical error is 


less than the size of the data points. 
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4.5 ea Cap) a05 compound nuclear reactions. 


Table 4.4 lists the energy levels of coe and 
is taken from En67. The levels are numbered according to 
the scheme of Be6é5. 

ag Gado) tloa reactions, many states in tloa 
are formed nearly exclusively by compound nuclear reactions. 
These states are characterized by low - cross sections and 
pelatively isotropic angular distributions. Cross sections 
for five such states (levels 7, 8, 9, 10 and 16) were 
measured between 50° and 170° in 10° increments, at energies 
from 4.50 MeV to 5.50 MeV in 2 keV Steps. As mentioned 
in 3.2, the cross sections were automatically averaged over 
90 keV energy increments, then averages over larger energy 
increments were done manually. The effect of the size of 
the energy averaging interval on the cross sections was 
imvectigated byjlaveraging Over 250) keV G10) 500) kev Cs 201) 
and 1000 keV (>40r). 

Figure 4.11 shows four angular distributions for 
state 9 obtained by averaging over different 250 keV inter- 
vals. Since the mean energy for each interval is different, 
absolute cross sections cannot be compared. However, we 
can compare the shapes of the distributions. It is fre- 
quently assumed that averaging cross sections over energy 
itcervale of 101 or 20% 268sutilelent to remove fluctuations, 
making the results directly comparable to Hauser-Feshbach 
theopy, Ve.c. st/0%, «booda). » bhegdistrricions an) figure sited 
are averaged over >101; therefore , if the above assumption 


ie correct, we would expect very little change anjshape for 


a” meh y cgi 8 


Reet) a ies ae 


pevadmun 27s iit ot vane 
aye az essete enem remoltosey 5th Lae ee 5 
,2holtaest teglsun Uatoamos yd cee earitickin ylasan bear? om | 
bas efoltoaa eso%o wol yd beattstosiBds ace astete eeodT - 


snoitses aBerD ,anotsiudidteth rsivans oigostesk ievenelaea 


_ 


Se 


evsw (df brs Of .% #8 .\ efevel) aspete douse ovit got 
asigisie ts, .etoasmyetons °Of ni °OVL bas 02 neewitsd bomuae 
bsnoisnsm| eA .eqeta Vet 8 At Ve" 02.2 or Vem 0e.# mort 
TSVo begevers va fecbpemotus biaw ‘erortbee: seus. StF Sek at 
YS IS9 regret YWsvo deyeisve nett ,etnomsront yaesae Yes) Oe 4 } 
to ssia-sdt to tosis sit ,vIfuvhem snob ouisw sinemegonin 7 
SBW enolrore 2201s edt no Lsevretnt anbasiwve yoIsne ens 
(10S<) Vox O62 ,(10Ls). Vet O28 tevo gninsisve Yd Setsgeteeee 
~CTOH#<) Vor O00L Brs- 
Yot enottudinterb talvens Yuet awode Lf. savgtd 


: 
i 7 
(3J 


-isini Vsxt O¢cf Fnomelith csyvo anissisve vd beatstdo ® — ; 


eraasStiscb et [feyrsial rossi rol yexene asem ant sont. -elsv : 

i a 
ow ,tsvewol .bsyaegnon ed tonns> anoitsse aaoto séploads a7 
“932 ef 31 Jenottydiutarb isis tc asqede sshd StBqMOD 159. 7 


YE9SNS isvo enoisose 22010 gatgsvovs vant Samuses vEsmOUP.. _ 


e@noivsutovl? svomsa.ot+ tneiottivne er 1NS 1o 10f Yo efavessni 
paadateteresisyat ot yitoeutb ativess sit gmk sa 
CN vaugit nb oathacds. at “thea 2ONI8! ..9) on 
mobigaveas svads ors" % = arenes glOLe Fave: begee ". 
ne Sgnde at onto 3 b #11 av. tons bulge ‘oom 


Table 4.4 


First 20 energy levels of *lcoa (En67) numbered according 
to scheme of Belote et al (Be65). 


Level number Excitation Energy 
. (MeV) 
0 0 
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Zz 2 WO 
3 2.463 
4 ee BIR 
2) 206 
6 22070 
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Figure 4,11 . 
*9ca(d,p)ticas: 9 angular distribution 


(state 9). Averaged data, averaging interval = 250 keV. 


Figure 4.12 
As in figure 4.11, but with averaging 


interval = 500 keV. 


In both cases, the statistical error is less 


than the size of the data points. 


Figure 4.13 

As in figure 4.11, but with averaging interval = 
1000 keV. The distribution is plotted twice; the upper 
curve with a scale comparable to that in figures 4.11 and 
4,12. Hauser-Feshbach predictions for several spins are 
shown. If no error bars are shown, the statistical 


error is less than the size of the data point. 
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Figure 4,14 
40 ee : ‘ ; 
Cad, pie =Ca*—angular-distrabutions-for state; / 
(A1lea2-88) , state 8 GEC 20 2), state 10 (ttea?**9) ana 
state 16 (41ca3-68), Averaging interval = 1 MeV, mean 
bombarding energy = 5.00 MeV. Hauser-Feshbach calculations 


are given for states 7, 8 and 16. 
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a different mean energy. This is obviously not the case. 
This suggests that the averaging interval may be too small. 
We would also expect symmetry of the angular distribution 
about 90° if the state is populated entirely by the compound 
nuclear mechanism. The observed asymmetry could arise from 
the energy averaging interval being too small, or from 
small contributions from the direct reaction mechanism, or 
Dove he 
Figurey12 shows the results for an averaging 
interval of 500 keV (220r). The shape still changes with 
energy and there is asymmetry about 90°. 
In figure 4.13 the data have been averaged over 
the entire 1 MeV range of the data. The cross sections are 
much more symmetric about 90°. Since the averaging interval 
is now 240r, it is suspected that fluctuations must have 
been averaged out, and that the remaining asymmetry is due 
to a small direct’reaction contribution to the. cross sections. 
A similar suggestion has been made in Wi66, when asymmetries 
were found in a state formed by the compound nuclear mechanism 
in the a Ol ones reaction, even after averaging over 50°. 
The other compound nuclear states considered (i.e. 
7,8, 20 andpl6)ibenavyedianga Similar manner togstate 9: 
the 250 keV and 500 keV averaging intervals seemed to be 
inadequate. The angular distributions of these states 
for a 1000 keV averaging interval are presented in figure 4.14. 
Hauser-Feshbach calculations were made for every 
state considered in eter whether the state was formed 


primarily by the direct or compound nuclear mechanism. 
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Recall that in equation 2.8 , the Hauser-Feshbach predic- 
tion for the differential cross section as a function of 
angle, there is a sum in the denominator of the transmission 
cocificients to all possible decay States. In this study. 
decays were possible leaving 1208 excited at up to nearly 
2 MeV. We venly have information on the levels cof tloa 
up to 6.8 MeV, and only a few spins and parities are known. 
In such cases, the sum over transmission coefficients to 
discreet states can be replaced by an integration, but for 
this we need to know the density of states as sae function 
of their spin and excitation energy. The usual method of 
determining the density of states is as follows: 
Ci) “the density op 7G) sot States, Ol spill Oiled 
nucleus is assumed to be factorable into a term de- 
pendent only on the excitation energy E*, and 
another term dependent only on the spin J. The 
usual form for this is 
ps(E*) = p(E*) (20+1)exp es 
where o is a parameter called the spin cut-off 
parameter. o was determined by fitting the above 
distribution to the distribution of levels of known 
spin in ae (Gimied 5 Uy EGY) 
(2) as mentioned in 5.1, p(E*) is determined for 
a particular nucleus from a plot of log NCE) vs YE, 
where NCE) is the number of states below an excita- 
tion energy E. Figure 4.16 is such a level density 
plot Lon 4lea, The lonepoint at 8.65 MeV is ob- 
tained from S-wave neutron resonance measurements 


on '[Ca, including the effects of our spin 


Ubuse aia nt 


qinssn ot qu Ts falas nat? sriveat efuieson oe 


55!" 30 efsvsl sit no nottemroin ays ylao ow Vem St " ‘ 
snwond $¢e;estviteq bre enkge wats ylae bn& , Vom, 8.8 oF qu - 
_— 


ot dtngickitesoo nokeaimensat tsVo mus Sit ,esebo dover 


yol stud ,mortsagstar ns vd Beosigqes od neo estate Jastoeib | 


soiftonut 8 e6 estazsa to yiiensb eft wordt oF been aw eltrit ; 


40 bortom [suey ed? | eygteie aolttatioxe bas nige shed? te 


:ewollot ea 2b estst2 to yireash oft gninimuszeb 


5 to Ll nige to estste to (*3)1q vitensp edt (1) 
-sb mies s otat siiisyotss? sd of bsmyees ar eysioun 
bob .*2 vetand noitsesisxes sit no ying tnsbnsq 
3iT i mige sdt mo yino tasbasqesb mist sedtons 
et eint-xoi mrot Ladev 

ses qxs(L4tS) (#d)q = (4D) z0 
Tio-tus nige sit bsiles wTemsrsq 6 wi o. s1sdw 
s¥Vods sdf sartsiit yd benimistebh asw o  .tsismaisg 
fttwonk Lo efevel to nottudiareis any ot noLtudratarh 
(2f.¥ .gtt) so4* nk mage 
tot beninvistsb ef (*4)q ,{.,2 nit bonottnem es (Ss) 


a ev (4) gol to tolq 6 moa? eusioun xelyottasq s 
-Btioxs a= wolsd estate to asdmun sz et ee sotodw 
ne mot 


etiensh fevel 5 dove eh Afi otwgit a 
wdo 2i VeM: 2956 to ii nor it i 


— 


Figure 4,15 
Spin distribution of states with known spin 


in 4loa, The theoretical curve is calculated from 
a (20+l)exp (-J0J*1) 
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distribution. 

Hauser-Feshbach calculations for the +Uca(d,p)ttca 
reaction have been done previously at this laboratory by Roy 
and Davison (Ro69), using the program HAUSER. They compared 
the results of these calculations to the experimental mea- 
Surements of Belote et al (Be65) at E, = 7.0 MeV. Many of 
the observed non-stripping states were adequately explained 
by the statistical model, and spin assignments were made. 
The level density of *lca was fixed very accurately by 
requiring the statistical model calculations to give the 
proper absolute magnitude for many non-stripping states. 

In this study, similar calculations were done for 
the 1 MeV averaged *0Ca(d,p) tea data. Several points con- 
cerning the calculations should be made: 

(1) The program HAUSER was used. 

(2) The calculations were insensitive to the 

deuteron optical model parameters used. Two sets 

of deuteron parameters were used: those obtained 
in this study from the 1 MeV averaged data, and 
those from Schwandt and Haeberli (Sc69). Both 
sets gave quite similar results. 

(3) The absorptive potential for the deuteron 

used in the Hauser-Feshbach calculation includes 

the effects of deuterons absorbed from the beam 
by direct reactions as well as those absorbed 

by compound nuclear reactions. This introduces 

an error into the predicted cross sections, 


since we are only considering compound nuclear 
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reactions in the calculation. The usual correc- 


tion is calculated in the following way: 


= [—tot="D}}q 


SCN cn? 


Here Soy is the cross section given by the com- 
pound nucleus program, ocy is the corrected value 


to be compared with experiment, Len tie scone. 


oo G 
reaction cross section, and op is the total direct 
interaction cross section. 

In this investigation, Stot Was calculated 
by the program SNOOPY, and oy was estimated by 
DWUCK calculations using the experimental spectro- 
scopic factor (see section 4). The (d,n) direct 
cross section was calculated from the data of 
Leighton et al (Le68). This gave a correction 
igs! Bamfone) one ne 
(Y) "The level density used for)! Ca waa that 
obtained by Roy and Davison in Ro69. Other level 
densities were obtained from plots of N(E) vs. 

/%, as explained previously. 

(5) Since energy averages of the experimental 
data were taken over considerable energy intervals, 
i.e. up to 1 MeV, HAUSER calculations were done 
over the same interval in 25 keV steps. The 
theoretical cross sections were then energy 


averaged before comparison with the energy 


averaged data. 
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(6) There has been no arbitrary normalization 
of the HAUSER predictions. Absolute cross sections 
are calculated, and are directly comparable to 
the data (after correcting for direct interactions). 
The results of these calculations are partially 
given in figures 4.13 and 4.14, where they are compared with 
states in 'lca that are thought to be formed primarily by 
the compound nuclear mechanism. As mentioned previously, 
despite averaging of the data over >40r, and careful 
attention to background estimation, electronic dead time 
and target non-uniformity, the experimental angular dis- 
tributions of these states are not symmetric about 90°. 
They are slightly peaked in the forward quadrant, suggesting 
a small admixture of the direct reaction mechanism. 
In four of the five states considered (7, 8, 9 and 
16) the angular distribution is fairly symmetric, and H-F 
calculations were done in these cases. Peak 10, however, 
is strongly peaked, and no calculation was attempted. 
State 7 has J” = (7/2*) (Jo69). HAUSER predictions 
fOPpuspins 7/2 and’ 9/2 both give reasonable mits to wthe data. 
State 8, measured as J" = (7/27) (J069), is fairly 
asymmetric about 90°. HAUSER calculations for spins 7/2, 
9/2 and 11/2 give reasonable values for the absolute cross 
sections, if not the shape of the angular distribution. 
Johnson (Jo63) found state 9 to have spin 3/2 or 
5/2. The 1 MeV averaged data for this state give an angular 
distribution that is fairly symmetric about 90°. HAUSER 


predictions for spin 3/2 and spin 5/2 both give acceptable 
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4.6 rea Ca oe +loa direct reactions. 


Cross sections for states in *1ca formed predominantly 
via the direct reaction mechanism were calculated theoreti- 
cally using the DWBA code DWUCKT. All calculations used a 
finite range parameter g = 0.60 fm-! and nonlocality ranges 
of 0.50 fm for the deuteron, and 0.85 fm for the neutron and 
the proton. The optical model parameters used in these cal- 
culations are given in table 4.5. The deuteron parameters 
are those obtained in this investigation from the elastic 
scsttering data and the proton parameters are those of 
Greenlees and Pyle (Gr66). For the neutron, the depth of 
the central well is adjusted by the program so that the 
binding energy of the well is equal to the experimental 
separation energy, the spin-orbit radius and diffuseness are 
thewsamesiass ifort ithes cent radawelden.s ¢bhus: mthens ,are sohree. unde— 
termined parameters: the radius rso, the diffuseness aso 
and the depth Vso of the spin-orbit potential. These were 


taken from the work of Perey and Buck (Pe62). 


Table 4.5 
Optical model parameters used in DWBA calculations. 


Units MeV or tm. 
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Sample calculations were also done using the 
deuteron parameters of Sc69. The results were not signifi- 
cantly different, and are not presented here. 

The states at10.0, 1.943) 23010 and 22463. Mev 
excitation (states 0, 1 and 3) were selected for detailed 
Consideration. These states have relatively high (d,p) 
cross sections and are well separated from adjacent peaks. 

In the shell model, we picture 49Ca as a very 
inert nucleus, with neutrons and protons filling all shell 
model levels up to the 1 d3/2 level inclusive. Ideally, 
we could think of 41ca as a single neutron orbiting in a 
potential created by the inert +[2ca core. The ground state 
of tlca would then be a neutron in the 1 £7/2 shell model 
orbit plus the inert core.- As explained in section 2.2, 
DWBA calculations are done assuming this simplified picture 
to be true. The ratio of experimental cross sections to 
DWBA predictions gives the spectroscopic factor, which is a 
measure of the validity of the above model, i.e. a measure 
of the single particle character of the state. Experimen- 
tally, weuxnow that it 16 veryodiftaicult to excitesprocons 
or neutrons in the +9ca nucleus (it is "doubly magic"), so 
it should be a very inert core in 4lCa. Thus, we expect 
that at least the ground state, and possibly some excited 
states, should correspond very well to the model of a 
single particle orbiting an inert core. Because of these 
properties, 4lca has been rather extensively investigated. 
In Le68, the 40ca(d,p)ttca reaction was done at 5.0 MeV, 


and angular distributions for states 0, 1 and 3 were 
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obtained, These distributions are given in figs. 4.17, 
4.18 and 4.19. DWBA calculations were performed, but fits 
were poor in the backward quadrant, particularly for ene 
ground state. Part of the purpose of the present investi- 
gation was to see if the poor fits at back angles (where the 
compound nucleus mechanism is most important in a stripping 
reaction) could be improved by obtaining energy averaged 
cross sections, thereby eliminating the effects of compound 
nuclear fluctuations. DWBA calculations were done from 4.5 
to 5.5 MeV in 25 keV steps, and the results were energy 
averaged. The angular distributions obtained in this way 
were slightly different in shape and magnitude from those 
calculated at the mean energy. 

Figures 4.17, 4.18 and 4.19 give the unaveraged 
data of Leighton (Le68), and the 1 MeV averaged data and 
DWUCK calculations from the present investigation. Compound 
nuclear contributions were calculated using HAUSER, and 
subtracted from both sets of experimental data. 

For the ground state (state 0), the 1 MeV average 
greatly affects the experimental cross sections at angles 
greater than 90°, and improves the agreement between theory 
and experiment. The resulting spectroscopic factor is 0.83. 

The averaged distributions for states 1 and 3 have 
different shapes from those of Leighton at back angles, but 
there is no improvement in the fit of the DWUCK calculations 
to the data. There is a difference in the absolute cross 
sections also. The spectroscopic factors obtained for 


states l and 3 in the present investigation are 0.37 and 0.13. 
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Figure 4.17 
*Qeacd.p) teal 4 angular distributions: 


-- 1 MeV averaged data from this study, E 5.00 MeV 
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Figure 4.18 
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Cala. oe) Ca angular distributions: 


-- 1 MeV averaged data from this study, E 5.00. MeV 


mean 


-- Data of Leighton, Eg = 5.00 MeV 


-- DWBA fit to data, S = 0.37 
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Figure 4,19 
oOGa(d.p) Cate angular distributions: 


-- 1 MeV averaged data from this study, E 5.00 Mev 


mean 


=-—Data Of Leighton, Eqg= 9.00 MeV 


-- DWBA fit to data, S = 0.13 
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Figure 4.20 

Yield curve of “Sea (d,p) leaving tlca in ground 
state, first and third excited states. Cross sections 
are averaged over 50 keV. Measurements were made at an 


angle near the stripping peak. 
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The values for the spectroscopic factors are 
different from those obtained by previous investigators. 
Table 4,6 summarizes the values obtained from this Scudye 
those of Leighton et al (Le68) obtained at 5.00 MeV, those 
of Lee, Schiffer et al (Le64) which are averaged over the 
Pange 7 "12 MeV, “and tnose Of Senet al (Sel 0) ebtaimed 
aie ew TeV. 

Table 4.6 


Spectroscopic factors for states in tloa 


otate 0.0 IWSSENS) 2040 
Source 

This study 0.83 ORS Ona 
Leighton et al 0.68 0.2.50 Rake) 
Lee et al 0.86 Olee/ 2 Owe 2 
Seth et al 0.7 6 Ojo Ones 


The work of Lee, Schiffer et al and Seth et al 
was done at high enough energies that the compound nuclear 
contribution to cross sections is negligible. We therefore 
believe the spectroscopic factors deduced in those studies 
to be better than those obtained at lower energies. 

Referring to table 4.6, we note that the spectro- 
secopie factors obtained in this study for the 1.95 MeV and 
2.46 MeV states in 4lca are in substantial disagreement with 
the values, obtained at higher wenergy. \Wejlurther nove sche 
disagreement between the spectroscopic factors obtained in 
the present study (averaged data, mean bombarding energy = 


5.00 MeV) and those obtained by Leighton (data taken at 
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One) eneroy-g.o. 0. Mev)s 


Upon inspection of figure 4.20, at least part 
of the reason for the discrepancies becomes apparent. 
Figure 4.20°15 a@ yield? cunve for States, Beandiettaken 
at angles near their respective stripping peaks. Since 
the compound nuclear contribution to the cross sections is 
relatively small (<300ub/sr) we expect the direct reaction 
mechanism to predominate. Thus, a smooth yield curve is 
expected. The yield for the ground state is indeed fairly 
Ssmootir-and *the “spectroscopic #actor "obtained. PeNclose to 
the high energy value. However, minor fluctuations are 
present. The other yield curves have much more structure, 
especially the one for the 1.95 MeV state. It was surprising 
to find such large fluctuations in the cross section of a 
direct state at the stripping peak. It implies that spectro- 
scopic factors deduced by measuring cross sections at one 
energy in this range must be viewed with some suspicion. 
In fact, if we accept the spectroscopic factors measured 
at higher energy as being correct, our method of measuring 
cross sections in’ small steps over a’range* of °) MeV and 
averaging them is not sufficient to remove the effects 
of the structure of the yield curve and give good spectro- 


scopic factors. 
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4.7 Summary and Conclusions. 

Gis) The average width of levels in 4256 at “15 MeV 
excitation has been found to Deis 720 keV. 

G2) Data for +[Ca(d,p)41ca reactions were obtained 
between 4.50 MeV and 5.50 MeV in 2 keV steps. Energy 
averaged cross sections were extracted from these data. 
Averaging intervals of >10r and >20r were found to be 
inadequate to average out compound nuclear fluctuations. 

An averaging interval of >40r still gave angular dis- 
tributions for states 7, 8, 9, 10 and 16 that were not 
symmetric about 90°. (These states are expected to be 
formed nearly entirely by the compound nuclear mechanism. ) 
The asymmetry was attributed to a small but significant 
contribution from the direct reaction mechanism, giving 

a peaking at forward angles. 

(3) Comparison of the 1 MeV averaged data (mean 
energy = 0.00) for states’ 0, 1 and’ 3 with unaveraged data 
obtained by Leighton at 5.00 MeV showed changes in relative 
and absolute cross sections. This Was NOU expec veduat 
forward angles where the direct mechanism predominates, and 
suggests that cross sections should be energy averaged 
before comparison with DWBA calculations. Comparison of 
spectroscopic factors obtained from this data with reliable 
spectroscopic factors obtained at higher energies showed that 
an averaging interval of 1 MeV is sufficient to remove the 
effect of these fluctuations on the spectroscopic factor 
obtained for the ground state of 4loa. The spectroscopic 


factors of the first and third excited states obtained from 
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the 1 MeV averaged data do not agree with the high energy 
spectroscopic factors. This suggests that an averaging 
interval of 1 MeV is not sufficient in these cases. 

In view of the unexpected large fluctuations in 
Ene Cross section of the first excited state at .its stripping 
peak, it is suggested that in the future, before accepting a 
spectroscopic factor measurement, yield curves should be 
taken on the stripping peak to establish that the yield 
varies smoothly in the energy range of interest. 
(4) The agreement between experiment and DWBA at back 
angles for the ground state of 4loa was improved if 1 MeV 
averaged data were used rather than unaveraged data. 
Gay) The above procedure did not improve the fit for 


the 1.943 and 2.463 states of ‘ca. 
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